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ABSTRACT 

We report the discovery of a transiting brown dwarf companion to KOI-205, a K0 main-sequence star, in a 11.720125-day period 
orbit. The transits were detected by the Kepler space telescope, and the reflex motion of the star was measured using radial veloc- 
ity observations obtained with the SOPHIE spectrograph. The atmospheric parameters of the host stars were determined from the 
analysis of high-resolution, high signal-to-noise ratio ESPaDOns spectra obtained for this purpose. Together with spectrophotometric 
measurements recovered from the literature, these spectra indicate that the star is a mildly metallic K0 dwarf with T eff 5237 ± 60 K. 
The mass of the companion is 39.9 ±1.0 Mj up and its radius is 0.81 ±0.02 Rj up , in agreement with current theoretical predictions. This 
is the first time a bona fide brown dwarf companion is detected in orbit around a star of this type. The formation and orbital evolution 
of brown dwarf companions is briefly discussed in the light of this new discovery. 
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1. Introduction 

The formation and evolution of brown dwarfs (BD) is not well 
understood. While it seems clear that these objects are too mas- 
sive to generally form by core-accretion in a circumstellar disk 
(Fig. 3 of Mordasini et al. 2009), a number of plausible for- 
mation theories have been advanced (see reviews by Whitworth 
|et al.|20O7t|Luhman|2012| l, but the role each of them plays in the 
actual formation process is not clear. The small number of BD 
companions to solar-type stars prevents reliable statistical stud- 
ies from being made, and although radial-velocity surveys have 
detected a few tens of objects with minimum masses in the BD 
regime (e.g. Sahlmann e t al.|20lT Diaz et al.|2012 1, the inclina- 
tion degeneracy does not allow us to be certain that these objects 
are actually sub-stellar. Since the mass function is expected to 
rise in the stellar regime, the inclination effect is more severe 
for BD companions ( Ho & Turner|20lT i. On the other hand, in 
recent years transiting surveys have provided a handful of bona 
fide BD companions (e.g. Bouchy et al.||2011a Anderson et al. 
[20TTT >. 

Most of the discovered short-period BD companions are or- 
biting around F-type stars. This led |Bouchy et al. ( 201 lb| l to ad- 
vance the hypothesis that BD companions could not survive in 
close orbit around stars with convective envelopes (T e g between 
5000 K and 6000 K), since the tides excited on them by the mas- 
sive companion and the progressive spin-down due to magnetic 
braking would lead to the reduction of the orbit size and to the 
eventual collision with the host star. Because F-type stars are less 



* Based on observations collected with the SOPHIE spectrograph 
on the 1.93-m telescope at Observatoire de Haute-Provence (CNRS), 
France, and with the ESPaDOnS spectrograph on the CFH telescope. 



affected by braking throughout their lives, Bouchy et al. ( 201 lb} 
suggested they might be the most common long-term hosts for 
massive companions. 

In this letter we report the discovery of the first short-period 
BD orbiting a K-type dwarf, KOI-205, and discuss its implica- 
tions for BD formation and evolution. This object was part of the 
SOPHIE follow-up of Kepler candidates described by |Santerne] 
etaT71 ( |20T2] >. Like other objects followed-up with SOPHIE, KOI- 
205 shows that a significant fraction of the Kepler Objects of 
Interest are not transiting planets. 



2. Description of the data 

2.1. Kepler photometry 

KOI-205 (KIC 7046804) was observed by Kepler from Quarters 
1 through 13, i.e., between May 2009 and June 2012, with a sam- 
pling of one point every 29.4 minutes (Long Cadence data, LC). 
Although Short Cadence data are also available for some quar- 
ters, only LC data were used to reduce computation time, and 
because the transit shape is sampled well enough in the phase- 
folded data (Fig.[T|. The light curve issued from the Photometric 
Analysis module of the Kepler pipeline was used. One-percent- 
deep transits are clearly visible by eye; they occur every 11.7 
days. The typical uncertainty in individual LC points is around 
230 ppm. 

For the transit modeling only fragments of the light curve 
around each transit were considered. Each fragment was normal- 
ized with a parabolic fit to the out-of-transit part, and a sigma- 
clipping at 3cr was performed to reject outliers. The contami- 
nation by nearby stars was corrected using the crowding values 
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Table 1. Target coordinates and apparent magnitudes. 



R.A. (J2000) 

Kepler mag" 
SDSS g' " 
SDSS r' " 
SDSS !" " 
SDSS z' " 



Kepler ID 
2MASS ID 
19 41 59.20 



7046804 

19415919+4232163 

Dec. (J2000) +42°32'16'.'4 



14.518 
15.19 ± 
14.47 ± 
14.23 ± 
14.12 ± 



0.04 
0.04 
0.04 
0.04 



2MASS-J * 
2MASS-H * 
2MASS-Ks b 
WISE-W1 ' 
WISE-W2 c 



13.095 ±0.021 
12.709 ± 0.024 
12.656 ± 0.021 
12.564 ± 0.023 
12.688 ± 0.026 



Note s. (n) From th e Kepler Input Catalogue. "'' ICutri et aL| <|2003) (c) see 
|Cutri et al.| ( [2012) 

from the Kepler archiv^] Due to the quarterly rotation of the 
spacecraft, the crowding value changes every three months, but 
recurs every four seasons. After contamination correction, it is 
obvious that transits from Season 2 are systematically deeper. 
We believe this is caused by an overestimation of the crowding 
factor, which is 12% in Season 2 and between 6 and 7 % in the 
other seasons. We therefore decided to analyze the light curves 
from quarters belonging to Season 2 separately, and fit a con- 
tamination factor independently. 

2.2. SOPHIE velocimetry 

Seven spectra of KOI-205 were obtained between February 2012 
and June 2012 with the SOPHIE spectrograph (Perruchot et al. 
2008, |Bouchy et al.||20T3] l. Observations were performed in 
high-efficiency mode, with resolving power A/AA ~ 40000, and 
exposure times ranging from 900 to 2700 seconds. The resulting 
signal-to-noise ratios (S/N) per pixel at 550 nm are between 7 
and 16. 

The spectra were reduced and extracted using the SOPHIE 
pipeline ( |Bouchy et al.|20 09 ), and the radial velocites were ob- 
tained from a Gaussian fit to the cross-correlation function with 
numerical masks corresponding to different spectral types. The 
charge transfer inefficiency (CTI) effect was corrected using the 
polynomial in Sant erne et al.| ( 2012| ). For faint targets such as 



these, the spectral orders at the edge of the wavelength range 
are usually too noisy, and adding them in the average cross- 
correlation function degrades the precision of the measurements. 
For KOI-205, 8 and 5 orders were discarded from the blue and 
red ends of the spectrum, respectively. The resulting radial ve- 
locities are plotted in Fig.fTjand listed in Table|2] 

No bisector effect (Quel oz et al.|2 001) or mask effect were 
detected in the SOPHIE data, supporting the sub-stellar compan- 
ion hypothesis. Bisector measurements are also given in Table|2] 

3. Host star 

The stellar atmospheric parameters were obtained using eight 
ESPaDOns spectra, acquired in service mode using the star + 
sky configuration, which provides a spectral resolution of 65,000 
on the night of June 30 2012 (program 12AF05). The resulting 
combined spectrum has a S/N ~ 90 at 550 nm per resolution 
element. The method is described in |Deleuil et al.| ( |2012| l and 
yields the values 7/ eff = 5210 ± 70 K, [Fe/H] = 0.27 ± 0.14 dex, 
and logg = 4.65 + 0.07 [cgs]. The projected rotational velocity 
was determined to be v sin i m — 2 + 1 km s _1 . A comparison with 



StarEvol evolution tracks (Lagard e et al.|2012] ) indicates that the 
star is a K0 dwarf, with a wide range of allowed ages, because of 
the slow evolution of low-mass stars. However, its slow rotation 
suggests that the star is older than the lower end of the range 
given in Table [3] 

4. Modeling of the data and parameter estimation 

The Kepler light curve and SOPHIE radial velocities were fit- 
ted together with the photometric measurements from Table [T] 
which provide a well-sampled spectral energy distribution 
(SED), to a model consisting of a star orbited by a dark compan- 
ion of a given mass and radius. The model and the Markov Chain 
Monte Carlo algorithm used to take samples from the Bayesian 
posterior are implemented in the PASTIS package, which is 
described in detail in Diaz et al. (2013, in prep.). Basically, 



we modeled the light curves with the EBOP code (EtzelJTWT 



Popper & Etzel 1981 1 using a quadratic limb-darkening law with 



coefficients interpolated from the tables of Claret & Bloemen 
(|201 !)• To deal with the distortion of the transit shape aris- 



1 http : / /archive . stsci . edu/kepler/kepler_f ov/search . 
php 



ing from the finite integration time (Kipping 2010) of LC data, 
we oversampled the model to five times the original sampling 
rate and re-binned to the LC sampling rate before comparing 
them to the data. The radial velocities were fitted to a Keplerian 
orbit. The model SED was obtained using an interpolation of 
the P HOENIX/BT-Settl synthetic spectral library ( jAllard et aT] 
2012) for a given T e s, [Fe/H], and logg of the host star. The in- 
terpolated spectrum is scaled to a given distance and corrected 
from interstellar extinction. The distance d and the color excess 
E(B - V) were included as free parameters in our model. 

The model therefore has 14 parameters, which are marked in 
Table [3] The mass ratio q has only a slight effect on the model 
and because of the normalization procedure of the light curve de- 
scribed above, it is not constrained by these data. It was neverthe- 
less included as a free parameter to take into account its effects 
on the error budget of the remaining parameters. Additionally, 
we included parameters describing the data: the out-of-transit 
flux of each of the two Kepler light curves, and a contamination 
factor for Season 2. Finally, the systematic errors in the radial 
velocities and in the light curves are modeled as an additional 
source of Gaussian noise with zero mean and variance o 2 .. The 
widths cry of each dataset were also fitted. 

Samples from the joint posterior distribution of the parame- 
ters were obtained using the Metropolis-Hastings algorithm (e.g. 
Tegmark et al. 2004) with an adaptive step size (Ford 2006), 
coupled with an adaptive principal component analysis to cor- 
rectly sample the posterior even in the presence of non-linear 
correlations. Non-informative priors -i.e., uniform or Jeffreys 
distributions- were used for all parameters except for T e ff, z, P, 
and T c , for which a normal distribution was deemed more appro- 
priate. For P, and T c the determination by Batalha et al. (2012) 
was used, but the width of the distribution was increased by an 
order of magnitude to ensure that it would not bias the results. 

Forty chains of 700,000 steps each were started at random 
points drawn from the joint prior, and 25 of them appeared to 
converge to the same solution. The remaining ones were stuck 
in regions of lower posterior probability, and were not consid- 
ered further. A modified version |Geweke| ([T992 ) diagnostic was 
used to determine when each chain had reached a stationary 
state. To ensure that the samples used to estimate the parameters 
and their uncertainties are independent, the chains were thinned 
using their correlation length (e.g. Teg mark et al.|2004 1 before 
merging. In total, over 20,000 independent samples of the pos- 
terior distribution were obtained in this way. The formal 68.3% 
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OGLE-TR-123b 



OpLE-TR-122 b 



T^ELT-l b 



CoRoT-15 b 



OiaQ5_b 



-L^C6343 C f 



Orbital Period 
< 5 days A 5 -10 days 



Mass [Mjup] 



Fig. 2. Mass radius diagram including transiting companions 
more massive than 10 M Jup . The dashed curves are the Baraffe 
et al. (2003 i isochrones for (from top to bottom) 0.5, 1, 5, and 
10 gigayears. The dotted lines are the isodensity curves for 10, 
25, 50, 75, 100, and 125 times the mean density of Jupiter. 
(Data from: |Siverd et al |2"012j 



201 l blal |Johnson et al [?2 01 1 



Triau d etal.12013 
Mlier et al. 2009 



2008H Winn et al .|2008||Pont et al.|2005||2006| ) 



Bouchy et al. 



Deleuil et al. 




Wavelength [A] 



Fig. 1. Data and best-fit model for KOI-205. For the SED, the 
best-fit spectrum is plotted as a solid black curve, and the inte- 
grated flux in each of the photometric bands is plotted as open 
circles. 



confidence intervals are listed in Tableland Fig.[T]presents the 
maximum posterior model. The median value for <jj was found 
to be roughly 100 ppm for Kepler light curve^] and around 40 
ms~' for SOPHIE velocities, possibly due to the low S/N of the 
spectra and to incomplete CTI correction. The light curve con- 
tamination of Season 2 is around 6%. 



5. Results and discussion 

Our analysis shows that KOI-205 has a BD companion with a 
mass of 40 Mj up and a radius of 0.8 Rj up , in a 11.7-day circu- 
lar orbit. The position of KOI-205 b in the mass-radius diagram 
(Fig. [2]) ag rees well with the theoretical isochrones by Baraffe 
et al.| ( 2053"] ) for a system with an age between 5 and 10 Gyr. If 



the age of the system were closer to 1 Gyr, which is permitted 
by the data, then the models would not reproduce the observed 
radius. However, given the slow rotation rate of KOI-205, it is 
likely that the age of the system is closer to 5 Gyr than to the 
lower end of the range reported in table [3] Indeed, an analysis 
of the pipeline-corrected PDC Kepler light curve yields a ro- 
tational period of around 41 days, in agreement with the spec- 
troscopic v sin i m value. Although KOI-205 b is the smallest BD 



detected as yet, its bulk density is similar to that of WASP-30 b 
( |Triaud et al.|2013[ [A nderson et al.|201 1) > and lower than that of 
LHC6343C (johnson et al.|201 1) . 

On the other hand, KOI-205 b is remarkable in two senses. 
First, it is the only object in the mass regime between massive 
BDs, such as WASP-30 b or CoRoT-15 b, and light BDs (or 
massive planets) such as KOI-423 b and CoRoT-3 b. Secondly, 
it is the only short-period (P < 10-15 days) transiting BD 
known to date in orbit around a K-type star. Indeed, the second- 
most massive object orbiting around a similar star is HAT-P-20 b 
( Bakos et al.|201 1| , with a mass of 7 Mj up , i.e., almost six times 
less massive. Moreover, radial velocity surveys have discovered 
only four non-transiting objects with minimum masses above 10 
Mj up orbiting stars less massive than 1 M in short-period orbits. 
Given that observational biases tend to favor the detection of this 
population, its apparent paucity raises questions about its origin. 

Formation of BDs can proceed by gravitational collapse of 
a molecular cloud in a similar fashion as stellar objects form. 
Recent simulations (Bate 2012) have shown that BDs can be 
formed in binaries with stellar primaries, and they reproduce 
quite naturally the paucity of BDs in close orbits. Another pos- 
sible formation mechanism is disk fragmentation. Numerical 
simulations (e.g. Thies et aL]|2010[ ) have shown that numerous 
BDs can be formed in the outer regions of massive circumstel- 
lar disks. The mass spectrum of objects formed around a 0.7 M Q 
star has a maximum around the mass of KOI-205 b (Fig. 4 of 
Stamatell os & Whitwort h 2009), so the measured mass is not 
only understood by but also expected from theoretical consider- 
ations. However, both processes form BDs no closer than 10 AU 
from the star, and therefore a mechanism is needed to migrate 
the newly-formed object to the inner system. 

Scattering by other objects formed around the same star (e.g. 
|Ford & Ra sio 2008) or the Kozai mechanism (Kozai |1962[ ) can 
bring the object to a closer orbit, possibly resulting in an ec- 
centric and non-aligned orbit. Subsequent tidal dissipation is ex- 
pected to have rapidly circularized the system (e.g. Matsumura 
et al. 2008). In these cases, the presence of at least one additional 



2 This value seems much higher than the typical value reported by 
|Gilliland et al.| pOl 1| >, but the noise metric they used greatly differs 
from ours. A rough transformation of cr, to their metric yields a value 
that agrees with the stellar noise distribution presented in that paper. 



object in the system would be expected, but it is also reasonable 
to assume that these additional bodies are either very far out or 
have been ejected. Unfortunately, the time span of our current 
RV data is not sufficient to probe these possibilities. 
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Table 2. Radial velocity measurements for KOI-205 



BJD 


RV 


O-RV 


BVS" 


bxp. time 


S/N/pix 


-2 450 000 


(kms- 1 ) 


(kms- 1 ) 


(kms" 1 ) 


(s) 


(550 nm) 


5984.6759 


12.144 


0.054 


-0.386 


1800 


10 


6015.6232 


18.724 


0.019 


-0.077 


2700 


16 


6064.4682 


16.055 


0.061 


0.193 


1200 


7 


6071.5153 


15.932 


0.030 


0.061 


1371 


11 


6100.5669 


14.414 


0.038 


-0.035 


1079 


9 


6103.3985 


11.213 


0.043 


-0.052 


900 


9 


6105.5132 


13.640 


0.036 


0.017 


1200 


11 



Notes. (fl) Bisector velocity span. 




Table 3. Orbital and physical parameters. 



n / n crit 

Fig. 3. Total angular momentum (top) and energy (bottom) of 
the system as a function of angular orbital frequency n = 2n/P, 
scaled to the critical values (see text). The solid curve is the coro- 
tation locus (n = Q s = 2n/P mt ) for zero eccentricity, and the 
dashed and dot-dashed curves are the n = 2Q S and n = 4Q, loci, 
respectively. The green dashed curve is the constant total angu- 
lar momentum locus. The position of the KOI-205 system as of 
today is marked by the square that is zoomed in the insets. 



The fact that the orbit is nearly circular suggests that the tidal 
interactions are strong in the system, especially considering the 
spectral type of the star and the relatively high mass ratio. In 
fact, as mentioned in the introduction, close massive compan- 
ions to late-type stars are not expected to survive engulfment. 
However, by taking the stellar spin rate from the spectroscopic 
v sin u and assuming that the rotation of the BD is synchronized 
with its orbital periocj^] it can be seen that the system is able to 
evolve toward a stable configuration (Fig. [3]l . Indeed, the total 
angular momentum of the system is higher than the critical value 
(Hut 1980) while the orbital frequency is higher than the stellar 
rotation frequency. Therefore, tidal evolution will lead KOI-205 
b to migrate toward the star but engulfment will be prevented 
because the system will reach a stable equilibrium state at coro- 
tation. This can be clearly seen from the bottom panel of Fig. [3] 
which shows an energy minimum of the constant angular mo- 
mentum curve at corotation. 

This picture is complicated by the magnetic braking of the 
star, which leads to a loss of total angular momentum. But even 
considering a Skumanich-type law for the loss of angular mo- 
mentum (Spada et al. 2011), and considering the extreme case 
where the escaped angular momentum is drawn entirely from the 
orbital angular momentum at constant stellar spin, the system is 



Period P* [days] 


11.7201248 ± 


2.1e-06 


Midtransit time T c * [BJD]" 


975.17325 ± 


1.2e-04 


Eccentricity e* 


< 0.03 




Argument of periastron co [deg] 


263 t 


61 
230 


Inclination f [deg] 


88.456 ± 


0.055 


Radial- velocity amplitude K* [km s -1 ] 


3.732 ± 


0.039 


Radius ratio k = R b /R s " 


0.09849 ± 


0.00049 


Mass ratio q = M b /M s " 


[0.0 - 


0.02f 


Center-of-mass velocity y* [km s -1 ] 


15.057 ± 


0.026 


Semi-major axis [AU] 


0.0987 ± 


0.0013 


Transit duration [hours] 


3.07 ± 


0.15 


Stellar density p s * [p ] 


1.550 ± 


0.073 


Effective temperature T cff * [K] 


5237 + 


60 


Metallicity [Fe/H]* [dex] 


0.14 ± 


0.12 


Distance d* [pc] 


585 + 


16 


E(B-V)* [mag] 


0.040 ± 


0.023 


Rotational velocity v sin i, [kms" 1 ] 


2.0 ± 


1.0 


Rotation period [d] 


40.99 ± 


0.50 


Stellar mass M s [M Q ] 


0.925 ± 


0.033 


Stellar radius R s [Rq] 


0.841 ± 


0.020 


Age [10 9 yr] 


[0.4 - 


8.3] rf 


Companion mass M b [Mj up ] 


39.9 ± 


1.0 


Companion radius R b [Rj up ] 


0.807 ± 


0.022 


Companion density p b [pj llp ] 


75.6 ±5.2 


Equilibrium temperature T cq [K] 


737 ± 


31 



Notes. w Fitted parameters in the MCMC algorithm <fl) BJDJJTC - 
2,454,000 < * ) 99% upper limit (r) uniform distribution; not constrained. 
(d) 99% confidence interval 



3 The characteristic synchronization timescale is ~ 500 Myr, which 
is short compared to the age of the system. 



expected to be stable on the scale of thousands of gigayears (Hut 
1980). Furthermore, if KOI-205 were a more massive G-type 
star with a five-day rotational period, the system would still be 
stable over more than 10 Gyr. 

This shows that stars with convective envelopes are capable 
of harboring massive companions in close orbits. Therefore, the 
lack of detections around late-type stars cannot be explained by 
their being doomed to fall into their stars. On the other hand, 
an inefficient formation process, selection biases, or insufficient 
statistics might be invoked. Whichever is the case, KOI-205 b 
will prove to be of great importance in constraining formation 
and evolution theories of brown dwarfs and massive extrasolar 
planets. 
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